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Abstract 

A new formula of experimental evaluation of the stress 
intensity factor was proposed depending on the area of the 
caustic. From series of experiments, which were performed 
in single edge-cracked specimens made of Lexan (PCBA), it 
was found that the values of the stress intensity factor Ki„ 
which were based on the diameters of the caustics, were 
varied with the load, the crack length and the caustic shape 
because the caustic become oval in higher loads or higher 
crack lengths. This means that the stress intensity values 
which are evaluated by the two diameters of the caustics (Dt , 
Di), are significantly differed. The accurate stress intensity 
factor Ki can be experimentally evaluated by the area of the 
formed caustic. According to the new formula of the stress 
intensity factor of the stress optical constants were evaluated. 
The present work illuminates the complicated problem of 
the stress intensity factors evaluation by the caustic 
diameters and establishes the new accurate of the caustic 
area approach for the evaluation of the stress intensity 
factors at crack tip and the stress optical constants of the 
transparent materials. 
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Introduction 

The optical method of caustics, introduced by Manogg 
(1964) and Theocaris (1970), has been shown to be very 
effective for the determination of the characteristic 
parameters of singular elastic fields, where the 
corresponding stress parameters in these particular 
regions are governed by singularities which render the 
solution of the difficult problem, if not impossible, by 
conventional stress analysis methods. The difficulty 
arises mainly from the fact that the highly strained 
region at the singularity is very small and therefore 
the information gathered by classical experimental 
methods is rather vague and inaccurate. 

According to the method of reflected caustics, a light 
beam impinges in the immediate vicinity of the 
singularity and is reflected from the front and rear 
faces of the plate (for the case of transparent materials). 


The reflected rays are collected along a singular 
surface, which is strongly illuminated. A reference 
screen, placed some distance from the plate, intersects 
this surface and yields a singular curve, the caustic, 
which is, for all cases studied up to now, a generalized 
epicycloid. In this way, and by applying simple laws 
of geometric optics, that is the reflection laws, the 
singular stress field is transformed into an optical 
singularity, the caustic, which yields all the necessary 
information for evaluation of the stress singularity. 

The method of reflected caustics was used for the 
solution of several elasticity problems of particular 
interest in engineering applications. The method was 
applied primarily to cases with singularities, such as 
those appearing in cracked plates under plane stress 
conditions; to contact problems; and to problems of 
multiwedges (composite materials). However, the 
method works equally satisfactorily in elastic 
problems with any type of stress concentration, not 
necessarily including singularities. In such cases the 
caustic is generated from a deformed boundary 
instead of from a singular curve in the interior of the 
plate (initial curve). For a review of all these 
applications the reader is referred to a review paper by 
Theocaris (1974, 1977, 1978), Rosakis (1983), Kalthoff 
(1987), Papadopoulos (1993, 2005, 2011). 

The method of caustics was theoretically developed 
for elastic homogeneous and isotropic materials. For 
virtual materials, the caustic shape dependes on the 
material behavior at crack tip at the board of which is 
defined by the initial curve. Also, the caustic shape 
dependes on the plastic zone, core region, which is 
developed at crack tip because of Poisson's 
phenomenon. The shape of the caustic can be changed 
by the specimen surface curvature (Theocaris 1976), 
(Tomlinson and Patterson 1999). 

The experimental method of caustics was applied to 
evaluate the stress intensity factors by the caustic 
diameters. The values of the stress intensity factors, 
which were calculated by the diameters of the caustic, 
( Z), transverse and D, longitudinal diameter), which 
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its ratio is D t / Dj =1.056, are equal, K^Dj) = Kj(D t ) . 
While, in higher loads the ratio of diameters is 
D ,ID, <1.056, which means D, < Dj , the calculated 
values of the stress intensity factors are 
unequal, K I (D t )> Kj(D t ) . Then, a new formula for the 
evaluation of the stress intensity factor which was 
based on the caustic area, was applied, Papadopoulos 
(2011). By this formula, a new approach to evaluate 
the stress optical constants, Badalouka et.al (2011). 
Recent work, Unger et. al. (2005) and Hedans et. al. 
(2008), has demonstrated that the method of caustics 
can be used to study the plasticity around the crack tip. 

The caustics dependes on the stress optical constants. 
The stress optical constants of PCBA (Lexan) and 
PMMA (Plexiglas) were calculated by the strain gauge 
and interferometric techniques, Raftopoulos et. al. 
(1976). In a series of recent publication, Konsta- 
Gdoutos (1996 a.b. 1999) was reported resulting from 
the character and the extent of a three-dimensional 
region surrounding the crack tip by the method of 
caustics. According to these publications the stress 
intensity factors can be determined even for the case 
when the initial curve of the caustic lies in the region 
where the state of stress in three-dimensional cannot 
be approximated as plane stress or plane strain. Also, 
in these publications nomograms of stress optical 
constants versus an empirical factor of triaxiality were 
given. In this work the stress optical constants were 
experimentally determined by the reflected caustics. 
The triaxiality factor can be experimentally 
determined by the stress optical constants variation 
according to nomograms of Figs 7 and 8 of Ref. 
Konsta-Gdoutos (1996b). 

Method of Reflected Caustics 

The optical method of caustics is able to transform the 
stress singularity into an optical singularity, using the 
reflection laws of geometric. For divergent light beam 
the reflected light rays from front, (f) and rear, (r), face 
of the plate form two caustics, the caustic (f) and the 
caustic (r). For a cracked isotropic elastic specimen, the 


parametric equations of the two 

caustics are 

Papadopoulos (1993): 


X r,f=A m r o(* os6»±G C os^) 

(1) 

Y rj = V 0 ( sin6,± ! sin y) 

(2) 



with 


C rJ = 


£z a dc r jK I 

A. m y[2x 


(3) 


with: 


Z + Z- 

1 = 


(4) 


where z o is the distance between specimen and 
reference plane, d is the thickness of the specimen, 
is the magnification ratio of the experimental set- 
up, Z ; - is the distance between specimen and light 
beam focus, s = 2 for caustic (r) and s = 1 for caustic 
(f) and Kj is the stress intensity factor for the mode-I 

stress state. The c r and c r are the stress-optical 
constants of the material. The stress-optical constant 
Cj is given by: 



(5) 


where V is the Poisson's ratio and E is the modulus 


of elasticity of the material. 

The experimental stress intensity factor Kj is 
estimated by the relation, Papadopoulos (1993): 

2-j2n 


K,= 


3 sz od^! 2 c rJ 


/ \ 5/2 

Rll 

5,j 


= 1.67C 


/ \ 5/2 

' n \ 

u t,i 


S t,i 


( 6 ) 


with C = 


£z 0 dA^ 2 Crj 


where D, / are the maximum diameters of the caustics 
( D t transverse and Dj longitudinal) and S t j are the 
correction factors with: 


S t =3.1702 andS l =3.00 (7) 

If the caustic is about circle. Fig. 1 (theoretical castic), 
the diameters ratio is: 

— = — = 1.056 (8) 

D, 5, 

Then, the K I values are: K l(l) ( = hC l(l)) . If the caustic 

becomes oval, Fig.l (experimental caustic), ^1.056, 

the Kj values are difference , K l(l) } ^ K I(D ) . In this 
state, the relation between the values of the K r is: 


where r is the radius of initial curve of the caustics. 
This radius is given by: 


K 
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Fig.l illustrates the theoretical and the experimental 
caustics at notch tip for notch length a = 0.069/7? for 
the stress a = 2.42 MPa . 



FIG. 1 THEORETICAL AND EXPERIMENTAL CAUSTICS IN 
LEXAN (BCBA) SPECIMEN WITH STRESS a = 2.42 MPa , FOR 
NOTCH LENGTH a = 0.069/77 

Because of the shape of caustic is about ellipse, its area 
is given by the relation: 

A = ^D t D, (10) 

According to caustic theory, Papadopoulos (1993), the 
diameters of the caustic are: 

£>,=3.1702Vo. A=3.00Vo (U) 

Then, from relations (3), (6), (10) and (11), the stress 

intensity factor Kj becomes: 

K, = 0.1358CA 5/4 (12) 


the caustic, lies between the two others. From these 
curves a linear variation and values coincidence of the 
stress intensity factor up to cr = \ 2MPa can be 
observed. After this stress, a non linear variation and 
values deviation can be observed. This means that the 
stress intensity factors values in this stress region, 
which were calculated with the area of the caustic, are 
more accurate than the others. 



FIG.2 VARIATION OF THE STRESS INTENSITY FACTORS 
(NORMALIZED TO C) CALCULATED BY THE DIAMETERS, 
K I(D }j Kj (D ) AND AREA, A , OF THE CAUSTIC vs THE 

TENSILE STRESS FOR THE CRACK LENGTH a = 0.040/77 

Evaluation of the stress optical constants 

The c r and c j are the stress-optical constants of the 
material. The stress-optical constant c, is given 
by: c. f -v / E , where v is the Poisson's ratio and E is 

the modulus of elasticity of the material. The stress 
optical constant c f can be calculated according to 

method of ref. Theocaris (1979) or by the method of 
strain gauge. The plotting of the two caustics for 
divergent light beam is illustrated in Fig. 3. 

From equations (1) and (2), according to Fig. 3 for 
0 = 0 , we get: 


Then, the energy release rate or the J integral is given 
by the relation: 


Gi 



0.0184C 2 d 5/2 


for plane stress state 


(13) 


Fig. 2 illustrates the experimental stress intensity 
factors (normalized to C) calculated by the diameters, 
^HD),^i(D,) an d area. A, of the caustic versus the 

tensile stress for the same crack length a = 0.040/7/ . 
From this figure, it appears that the curve of the stress 
intensity factor, which was calculated by the area of 


c, =- 


2[5-8.18Ajc/d ( ^ 2 ] 5/2 


(14) 


where Ar is the distance (BA) and A (r) is the area of 
the (r) caustic. By the relation (3) the ratio of the 
stress-optical constants c r ,Cf is estimated. Analogous 

relation to Rel.(14) from Fig. 4 can be written by 
combination of caustics (r) and (f) for convergent light 
beam. 


c^^-STSAx/^f 


(15) 
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FIG. 3 GEOMETRY OF CAUSTICS (r) AND (f), FOR DIVERGENT 
LIGHT BEAM 


Conclusions 

The proposed new formula, relation (12), which 
depends on the area of the caustic, was given more 
accurate and stable values of the stress intensity 
factors, as it was concluded from Fig. 2. The stress 
intensity factor becomes independent from the 
diameters of the caustic, which were in difficultly 
calculated if the caustic is not symmetric. 

The accurate values of stress optical constants for 
transparent materials can be experimentally 
determined by the reflected caustics. The stress optical 
constants dependes on triaxiality at crack tip region. 
From the variation of the values of stress optical 
constants and the nomograms of ref., Konsta-Gdoutos 
(1996), the triaxiality factor k can be experimentally 
determined. The stress intensity factors can be 
accurately determined from accurate values of the 
stress optical constants. 



FIG. 4 GEOMETRY OF CAUSTICS (r) AND (f), FOR 
CONVERGENT LIGHT BEAM 

The experimental stress intensity factor Kj is 
estimated by the relation (6). 

From relations (6) and (14), the accurate stress 
intensity factor can be estimated from Poisson's ratio 
v and the modulus of elasticity E or the stress-optical 
constant c, and the caustic area: 


K, 


0.1358 

7 o 3/2 

z 0 dA m C j- 
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